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Abstract 

We argue that the square of the norm of the Hubbard wave func- 
tion is proportional to the determinant of a matrix, which is ob- 
tained by linearization of the Lieb-Wu equations around a solu- 
tion. This means that in the vicinity of a solution the Lieb-Wu 
equations are non-degenerate, if the corresponding wave function 
is non-zero. We further derive an action that generates the Lieb- 
Wu equations and express our determinant formula for the square 
of the norm in terms of the Hessian determinant of this action. 
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Introduction 

Interacting many body systems share universal features. The best known 
example are the critical exponents which describe the thermodynamic prop- 
erties near critical points. Two systems that have the same critical exponents 
belong to the same universality class. In one space dimension there are many 
systems which are exactly solvable by Bethe ansatz [T, 0, Each of these 
systems represents a universality class. The Bethe ansatz provides a unique 
possibility to obtain exact results for these systems, which are not accessible 
by perturbative methods. 

In this article we address the problem of the calculation of the norm of 
the Bethe ansatz wave functions of the one-dimensional Hubbard model. 

In all cases known so far the norms of Bethe ansatz wave functions are 
proportional to the determinant of a matrix, which is obtained by lineariza- 
tion of the Bethe ansatz equations around a solution. This assures the non- 
degeneracy of the Bethe ansatz equations in the vicinity of a solution, if the 
corresponding wave function does not vanish, and thus provides an important 
piece of information about the mathematical structure of the Bethe ansatz 
equations. 

The first time a determinant formula for the norm of a Bethe ansatz 
wave function appeared in the literature was in Gaudin's article [ffl on what 
is nowadays called the Gaudin model. Gaudin's conjecture was later gener- 
alized to the XXX and XXZ spin-| chains |J and to the Bose gas with delta 
interaction [lj]. Due to the complicated nature of the Bethe ansatz wave 
function a proof was not available until the advent of the quantum inverse 
scattering method, which provided a powerful algebraic reformulation of the 
Bethe ansatz. 

The first proof of norm formulae for a number of important Bethe ansatz 
solvable systems was presented by one of the authors || in 1982. The proof 
given in M is rather general. It relies on the structure of the i?-matrix and 
the existence of an algebraic Bethe ansatz. It applies to the XXX and XXZ 
spin-| chains, to the Bose gas with delta interaction and to a number of other 
interesting systems. As a corollary follows the norm formula for the Gaudin 
model (for a direct alternative proof of the norm formula for the Gaudin 
model see 

In the following years the results of || were generalized and confirmed by 
different methods. Reshetikhin || generalized the proof of |J to the case of 
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a gls invariant i?-matrix, i.e. to systems with nested Bethe ansatz. Slavnov 
proved a formula for the scalar product of a Bethe vector with an arbitrary 
vector || . This formula is valid for models having the same R- matrix as the 
XXX and XXZ spin-| chains. It was recently confirmed by Kitanine, Maillet 



and Terras who used an algebraic approach based on Drinfel'd twists [10 



Proofs of norm formulae which do not rely on the algebraic Bethe ansatz but 
on the Knizhnik-Zamolodchikov equation or its quantized version were given 



by Reshetikhin and Varchenko JTT|] and by Tarasov and Varchenko fL2[ . In 
the latter paper norm formulae for systems with gl n invariant i?-matrices 
were obtained. 



Bethe ansatz solution of the Hubbard model 

The Hamiltonian of the one- dimensional Hubbard model on a periodic L-site 
chain may be written as 

L L 

^ = -EE (sVma + 4+i,a c i,*) + u E(«it - - \) ■ (!) 

i=i a=U i=l 

Cj A and Cj t a are creation and annihilation operators of electrons of spin a 
at site j (electrons in Wannier states), n^ a = c^ a c^ a is the corresponding 
particle number operator, and U is the coupling constant. Periodicity is 
guaranteed by setting CL+i, a = ci, a - 

The eigenvalue problem for the Hubbard Hamiltonian ([]]) was solved by 



Lieb and Wu [13] by means of the nested Bethe ansatz (for reviews see 
|3], [HJ). The Hubbard Hamiltonian conserves the number of electrons iV 
and the number of down spins M. The corresponding Schrodinger equation 
can thus be solved for fixed iV and M. Since the Hamiltonian is invariant 
under a particle-hole transformation and under reversal of all spins, one may 



set 2M < iV < L IJ 



We shall denote the coordinates and spins of the electrons by Xj and ctj, 
respectively, I, J.. Then the eigenstates of the Hubbard 

Hamiltonian (|IJ) may be represented as 



\N, M) 



L 

E E ^(^i,--- ) ^Jv;ai,...,a JV )c+ Jvajv ...c+ ai |0) .(2) 



JV • xi,...,x N =l ai,...,a]v=T,| 
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ip(xi, . . . , xn', di, ■ ■ ■ , cln) is the Bethe ansatz wave function. It depends on 
the relative ordering of the coordinates Xj. Any ordering is related to a 
permutation Q of the numbers 1, . . . , N through the inequality 

1 < XQl < Xq 2 < ■ ■ ■ < X QN < L . (3) 

The set of all permutations of N distinct numbers constitutes the symmetric 
group Sn- The inequality (^) divides the configuration space of N electrons 
into N\ sectors, which can be labeled by the permutations Q. The Bethe 
ansatz wave functions in the sector Q are given as 

i 4>(x 1 ,...,x N ;a 1 ,...,a N ) = 

^2 sign(PQ) <p P (a Q1 , . . . , a QN ) exp i k Pj x Qj ■ ( 4 ) 

P€S N V J = 1 / 

The function sign(Q) is the sign function (or parity) on the symmetric group, 
which is —1 for odd permutations and +1 for even permutations. 

Explicit expressions for the spin dependent amplitudes ^(agx, . . . , aQjv) 
can be found in [|TI|. They are of the form of the Bethe ansatz wave functions 
of an inhomogeneous XXX spin chain, 

M 

(p P (a Ql , . . . , a QN ) = Y A(X nl ,...,X 7rM )Y[Fp(X n i;y l ) . (5) 

Here Fp(X;y) is defined as 

\ iU/2 £t A - sin k Pj - iU/4 

A — sin kp y + lU /4 J=1 A — sin kpj + lU /4 

and the amplitudes A(X\, . . . , Xm) are given by 

A(Xu...,Xm)= n A "-T A I", iC//2 • (?) 

l<m<n<M Am An 

In the above equations yj denotes the position of the jth down spin in the 
sequence oqi, . . . , aQN- The y's are thus 'coordinates of down spins on elec- 
trons'. If the number of down spins in the sequence clqi, . . . , Oqjv is differ- 
ent from M, the amplitude <fp((iQi, ■ ■ ■ ,o>qn) vanishes. An expression for 
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<Pp(o,qi, . . . ,qqn) using the terminology of the algebraic Bethe ansatz can be 



found in [2C 



The wave functions (f|) depend on two sets of quantum numbers {kj \ j = 
1, . . . , N} and {A; | / = 1, . . . , M}. These quantum numbers may in general 
be complex. The kj and A; are called charge momenta and spin rapidities, 
respectively. The charge momenta and spin rapidities satisfy the Lieb-Wu 
equations (or 'periodic boundary conditions') 

= n "'" 5 '^:'^ , i=i,...,jv , ( 8) 

f-Jl Xi — sin kj + iU /4 
" A, -sin ^-1^/4 ji A; - A m - iU/2 

j^Xt-smkj + iU/A U Xl -\ m + iU/2 ' • [J) 

The states (§) are joint eigenstates of the Hubbard Hamiltonian ([!]) and 
the momentum operator (cf. fl?]]) with eigenvalues 



E = -2j2coskj + -{L-2N) , P = £ kj mod 2vr . (10) 
i=i 4 \i=i / 

The eigenfunctions (^) are antisymmetric with respect to interchange of 
any two charge momenta kj and symmetric with respect to interchange of any 
two spin rapidities A/. They are antisymmetric with respect to simultaneous 
exchange of spin and space coordinates of two electrons, and hence respect 
the Pauli principle. 

An action for the Lieb-Wu equations 

Let us rewrite the Lieb-Wu equations (§), (^) in logarithmic form, 

,ti liCZ-^-sm^j; ^ ^i(/-2(A,-A m ); ' ' 

Here the logarithm is defined in the cut complex plane, where the cut is along 
the real axis from — oo to zero. n C j in equation (|TTD is integer, if M is even 
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and half odd integer, if M is odd. Similarly, nf in fll2"|) is integer, if iV — M 
is odd, half odd integer, if N — M is even. 

In order to formulate our norm conjecture we shall introduce certain func- 
tions connected to the logarithmic form (|TTD, ( |l"2"D of the Lieb-Wu equations. 
We shall start with the definition 

e„ W = i/>i„(lf±!) • (13) 

In terms of this function the Lieb-Wu equations (|ll|), (0) read 

M 

k 3 L - %(h ~ sin kj) - 2-nn] = , (14) 
1=1 

N M 

Y^Q'uiXi-smkj)- £ Q' 2U (\i - \ m ) - 2Tmt = . (15) 

j=l m=l 

Here the primes denote derivatives with respect to the argument. The left 
hand side of these equations can be easily integrated with respect to the 
variables sin kj and A;, yielding the action 

TV 

S = sin kj + cos kj)L 

3=1 

N M ^ M 

+ J2J2®u(Xi - sin kj) - - Y @2C/(A«-A m ) 

j=l 1=1 1 l,m=l 

N M 

-27r^n J c sinfc j -2nY J n i X i ■ ( 16 ) 
j=i 1=1 

Thus, introducing the abbreviation Sj = sin kj, we can write the Lieb-Wu 
equations as extremum condition for S, 

dS dS 

— = , j = l,...,N , ^ = , / = 1,...,M . (17) 

We shall use the action S below in order to formulate our norm conjecture. 
A similar action was first introduced by Yang and Yang in the context of the 
Bose gas with delta interaction flIB|. 
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There is an interesting alternative way to write the Lieb-Wu equations. 
Let us define 



1 



M 



Xj = kj- -J^e^-sinfy) 



Mir 



1=1 



1 N i _ 

fl = M X 6 ^ A ' ~ Sil1 ~ ^? X] 2£/( A i - A r 



Af 



iV 



AT 



m=l 



(JV- M +1)tt 



(18) 



,(19) 



where j = 1, . . . ,N and / = 1, . . . , M. In terms of these new variables the 
Lieb-Wu equations (§), (!) become 



;^ L = 1 , j = l,...,N 



e mN = 1 , I = 1, . . . , M . (20) 



This suggests to interpret the Xj an d <^z as the momenta of charge and spin 
degrees of freedom. 



The norm formula 

The square of the norm of the wave function (|4]) is by definition 

L 

llVi 2 = X X! |^(xi,...,x iV ;ai,...,a i v)| 2 

xi,...,x N =l ai,...,ajv=T,| 

Note that the normalization in (||]) is such that 

(N, M\N, M) = I ip || 2 . 



(21) 



(22) 



We suggest the following formula for the square of the norm of the Hub- 
bard wave function (W), 



V Z 7 j=l l<j<k<M \ 

( d 2 S d 2 S \ 



M N 



u 2 



4(A i - A fc ) ! 



•det 



9s 2 <9s<9A 
<9 2 S <9 2 <S 



(23) 



V dXds d\ 2 I 
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Here M' is the number of pairs of complex conjugated k^s in a given solution 
of the Lieb-Wu equations. The determinant on the right hand side of ([23|) is 
the determinant of an (N + M) x (N + M)-matrix. This matrix consists of 
four blocks with matrix elements 

'd 2 S\ d 2 S r f L ™ U/2 } 




* 2 / ran 9s m ds n ' { COS k n ^ (£//4) 2 + (A/ - S n ) 2 J 

..,iV , (24) 



/d 2 S\ p 2 5 \ <9 2 <S U/2_ 

[dXds ~ [dsdX nm ~ dX m ds n ~ (<7/4) 2 + (A„ 



Sn) 2 



m = l,...,M , n = l,...,N , (25) 



'<9 2 S\ <9 2 S 



d\ 2 J dX m dX n 



■ f U/2 " U 



5 (C/4) 2 + (A„ - Sj y fri (U/2)' + (A„ - \,y 

V/2)»+L-A.)' ' '"•"= 1 --" M • < 26 > 

The norm is thus proportional to the Hessian determinant of the action 
S regarded as a function of the charge momenta kj and spin rapidities Aj. 
Recalling the formulation ([□]) of the Lieb-Wu equations we see that a solution 
is non-degenerate (locally unique for fixed values of n c - and nf), if the norm 
of the corresponding wave function does not vanish. In other words, all 
eigenstates of the Hubbard Hamiltonian (|l|) correspond to non-degenerate 
solutions of the Lieb-Wu equations. 

Another interesting form of the norm formula is obtained by express- 
ing the Hessian determinant in equation (|23D in terms of the momenta of 
elementary charge and spin excitations Xj an d <fi, 

I^J 2 = (-\) M 'N\L N N M 



(27) 




d(ki, . . . , kpj] Ai, . . . , Am) 
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The norm is proportional to the Jacobian of the transformation from the set 
of charge momenta and spin rapidities kj, A/ to the set of momenta of charge 
and spin degrees of freedom Xj, fi- 

Let us list our arguments in support of (|2~3|): 



Our experience with other Bethe ansatz solvable models || 0, |[ [L2 



shows that norm formulae for Bethe ansatz wave functions are generi- 
cally of the form (p3|). 

(n) It is easy to see that (^) is valid for M = and arbitrary N. 



(iii) We verified fl2|) for JV = 2, JV = 3 and M = 1. The calculation 
is lengthy and involves highly non-trivial cancellations based on the 
Lieb-Wu equations. 



(iv) We verified (|23| ) for arbitrary N and M in the limit U — > oo. This limit 
requires rescaling of the spin rapidities, Aj = AjU/2. The remaining 
non-trivial factor in the expression for the norm reduces to the case of 
the XXX spin- 1 chain, and the result of || can be applied. 

(v) We calculated (for arbitrary iV and M) the leading order term of the 
norm in the large L limit. This term is proportional to L N and fixes 
the prefactor in (|23|). 



Finally we would like to point out two applications of the norm formula 



23). The string hypothesis for the Hubbard model [19| together with the 



SO (4) symmetry [2(| predicts the correct number of eigenstates of the Hub- 



bard Hamiltonian |2l|. This fact is sometimes called combinatorial com- 
pleteness. For some Bethe ansatz solvable models there exists an alternative 
proof of completeness which avoids the string hypothesis P2| . An important 
ingredient of this proof is a norm formula similar to (|23|) . 

Over the past two decades a programme was developed for calculating 
exactly temperature correlation functions of Bethe ansatz solvable models 
0. In the thermodynamic limit the correlation functions are represented 
as determinants of Fredholm integral operators of a certain form. The de- 
terminant representation is a powerful tool for studying the properties of 
the correlation functions analytically. It enabled, in particular, the exact 
calculation of long and short distance asymptotics for a number of Bethe 
ansatz solvable models. The calculation of the norm of the Bethe ansatz 
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wave function is an important step towards a determinant representation for 
correlation functions of the Hubbard model at finite coupling. For the special 
case of infinite repulsion a determinant representation was recently obtained 
in @. 
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